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I propose a definition for a “shocking coefficient” 5 intended to make determinations of the degree of wave¬ 
form shocking, and comparisons thereof, more quantitative. This means we can avoid having to make ad hoc 
judgements on the basis of the visual comparison of wave profiles. 


I. INTRODUCTION 


II. SHOCKING COEFFICIENT 


The formation of ’’shocks” on wave profiles - regions with 
a step-like or otherwise excessively steep gradient - has been 
an interesting area of nonlinear optics for some time in 12. 
In particular, recent interest has centred on carrier shocking 
iiniffla where the underlying wave oscillations self-steepen, 
rather than just the wave envelope 12. Recently it has also 
been shown to have potential for practical - or at least ex¬ 
perimental - relevance, based on recent predictions for and 
mid-infrared materials with realistic dispersions ilEl. Con¬ 
sequently, in this work I will consider only this nonlinear op¬ 
tical situation, although the definition of shocking coefficient 
I propose may well have utility in other fields. 

Here I take a purist point of view in that a shock occurs 
(only) when the pulse profile has evolved to have an infinite 
gradient; this is distinct from a more intuitive criteria involv¬ 
ing the presence of a steep edge in the wave profile with (pos¬ 
sibly) some fast trailing oscillations. However, by insisting on 
an “infinite gradient” criterion, the notion of how “shocked” a 
given waveform or pulse might be is not so obvious. After all, 
any finite gradient, however extreme, is a long way from being 
infinitely steep. Further, even without that absolutist position, 
how do we decide on - and most importantly quantify - what 
features of the waveform are or might be shock-like? 

One notable exception to the difficulty of discussing shocks 
occurs in numerical experiments: one can use the numerical 
difficulties that occur when the pulse evolution reaches the 
limits of computational resolution. This is in essence what 
is done by the LDD (local discontinuity detection) method, 
as already used in an optical context by Kinsler et al. Q- 
However, the LDD method, although useful in its idealised 
context, cannot inform us as to the degree of shocking present 
in a given waveform, especially if it is not approaching the 
limits of numerical accuracy. 

In Sec. I^I consider a way to calculate a “shocking co¬ 
efficient” that can be used as an indication as to degree of 
shocking present in a waveform, and a as a benchmark for 
comparisons. Following that, in Sec. 01 present the results 
of some simulations using this coefficient, and I conclude in 
Sec.HV] 
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Here I propose a relative measure of how shocked a given 
pulse is, by the simple expedient of regularizing the gradient 
with the inverse tangent function, which can be used to map 
very large (or even infinite) gradients into a fixed and finite 
range. I define the shocking coefficient, for a point on a wave 
A(f) of (relative) gradient g, as 

5=-arctan(|g| - 1), (1) 

n 

where the Ijn factor ensures that an infinitely steepening 
wave returns 5 —?► 1, and the g — 1 ensures that an unchanged 
wave returns S — 0. If the wave has become smoother, the S 
will become negative, down to a minimum of 5 = —0.5 when 
the wave is flat, i.e. |g| =0. In what follows I choose the ref¬ 
erence value to be the maximum value of the initial gradient, 
as evaluated over some entire initial waveform Ao(f); i.e. 


8{t) 


dA{t) 

dt 


Max 


dAo{t) 

dt 


-1 


( 2 ) 


Note that here I consider waves specified as functions of time 
t as they propagate forward in space nMni, however there 
is no impediment to instead choosing to analyse wave pro¬ 
files defined over space using this scheme; just replace t above 
with any or all of.r,y,z as appropriate. Further, although there 
might be cases where it would make sense to use a scaling 
that varies depending on the position of a point on the wave 
profile, I do not address that case in this work. 

To demonstrate the plausibility of the definition in eqn. Q, 
let us first consider a simple example. We start with a wave 
of sinusoidal profile, A(f) = Ai sin(a)f), which has maximum 
gradient (O and intensity Aj®. Now imagine a convenient 
nonlinear process converting this into a multifrequency wave 
A'(f) with a weighted spectrum of N harmonics, all perfectly 
in phase. Note that to get every harmonic appearing in an 
optical context we would need to consider a suitable second 
order nonlinearity II2II2, although it is unlikely to give rise 
to an equal weighting between the generated harmonics. The 
new wave A' that is to be tested for shock-like gradients is 
therefore 


N 

A'{t) = ^ A„,sin(m®f). 

m=l 


(3) 


Assuming an equal weighting with all A„, =Ai^, summing of 
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the harmonics shows us that this has an intensity 

N .2 

/= ^A>® = ^iV(iV+l). (4) 

m=l 

which means that the amplitude of each harmonic - by con¬ 
servation of energy - must be 

Ai 

An= ! - (5) 

^\N{N+l) 

It then follows that the maximum gradient is 

^Aig = ^ Afjsm{mcot) = . * - . —N{N 4-1) 

m=i 2 

( 6 ) 

= (0Ai^/N{N+l)/2. (7) 

Defining = yjN{N + l)/2, and noting that for large N, the 
limit is that > N, the shocking coefficient for our wave¬ 
form A' is 

2 

5 = — arctan — 1]. (8) 

As an example, for waves A' with N up to 3,5,7,9 max¬ 
imum harmonics (i.e. A* values of 2.45, 3.87, 5.29, 6.71) 
we then find that S becomes 0.62, 0.79, 0.85, 0.89 respec¬ 
tively. On the basis of this admittedly simplistic calculation, 
we might expect than more than just the first few harmon¬ 
ics were needed if near shocking behaviour is to be indicated. 
This will, however, depend on the threshold value of S we 
chose to use to indicate near shocks. 

Another case is when only odd harmonics are produced, 
which (e.g.) occurs for the third-order nonlinear (Kerr) 
case typically considered in nonlinear optics and shocking 
m. Again making the simplifying assumption of equally 
weighted harmonics, the same process as above can be fol¬ 
lowed. By noting that the sum of a sequence of M odd 
numbers is simply M^, we find that the equivalent A* for 
this case is just M, and we can reuse eqn. (|^ accordingly. 
Now, waves A' with maximum harmonic numbers of 3,5,7,9 
(i.e. A* values of 2, 3, 4, 5) give rise to S coefficients of 
0.50, 0.70, 0.80, 0.84 respectively. 

Numerically we can estimate the S achieved by various 
power law fall-offs in the harmonic spectrum under the same 
conservation assumption as above. The results can be seen in 
fig-0 One interesting case is the 7 = 4/3 fall-off considered 
by Boyd 01 [la. As should be expected, however, once the 
power law exponent goes beyond —2 the potential for near- 
shocking gradients is greatly reduced. 

Although these sample calculations above are indicators 
of how this shocking coefficient works, the main drawback 
is that in realistic cases the energy in the harmonics builds 
up gradually, and we should not necessarily expect the wave 
energy to be redistributed among harmonics in some simple 
way; not to mention the likelihood of complicated phase re¬ 
lationships being present. Accordingly, in the next section I 



FIG. 1: Shock coefficient S for a multi-harmonic spectrum with 
Am generated from a fundamental wave at (Bj, assum¬ 

ing energy conservation. The solid line is the result for when all har¬ 
monics are present, the dashed line is for when only odd harmonics 
are present. 

analyse some nonlinear optical simulations which provide less 
artificial testbed. Nevertheless, these considerations based 
on hamonic content may be useful in setting threshold val¬ 
ues of S that indicate nearly shocked waveforms. For exam¬ 
ple, total conversion of a wave into its third harmonic gives 
S ~ 0.70 - so should we consider waves of comparable S val¬ 
ues “shocked”? 


III. SIMULATIONS 

The simplest simulations we might apply the shocking co¬ 
efficient of eqn. Q to would be for a dispersionless medium. 
Accordingly, in fig. 0we show the results for such a 
medium, based on data from Kinsler a which assumed an 
initial CW profile. Fig. |^a) shows the wave profile evolu¬ 
tion as the wave steepens, and fig. |^b) show the matching 
increase in the peak shocking coefficient. We also see that the 
increase in shocking coefficient in one part of the wave profile 
is counterbalanced by a decrease in the other part. Further, the 
same type of progression can be seen for the more commonly 
considered medium on fig. |^a,b), albeit with a different 
wave profile evolution. 

Next I turn to a less idealized situation, and present the be¬ 
haviour of S as indicated by simulations that also include dis¬ 
persive effects. Here, the simulation parameters have been 
chosen primarily for illustrative purposes. The material re¬ 
sponse is defined by a centre frequency (Oq = 2.26946 x 
a refractive index n = 1.5, group velocity Vg and 
a dispersion D 2 which gives a basic quadratic form of k{co) 
which is : 

1 9 

k{co) =ko + Vg{co-(Oo) +-D2 {co-(Oo) , (9) 

although note that below I will use rescaled parameters Vg = 
Vg X 10^ and d 2 = D 2 X 10^^. However, this originally spec¬ 
ified material response undergoes filtering and processing to 
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FIG. 2: Shock progression in a dispersionless P) medium, with the 
simulation data used to create from fig. 1(b) of Kinsler O. (a) Wave 
profile, where the solid, dashed, dash-dotted, and dotted lines are for 
propagation distances closer and closer to the shocking distance, (b) 
Shocking coefficient for the same data, which increases as expected 
and also indicates the location of the most shocked part of the wave. 


guarantee causal behaviour in line with the Kramers Kronig 
relations M- 

On fig. I^we can see how the shock coefficient 5(f) changes 
for increasing nonlinear strength. The behaviour is straight¬ 
forward and intuitive. For simulations where the dispersion 
overwhelms the nonlinearity-induced tendency to shock, a 
maximum S is reached after which the wave smooths out, and 
over longer periods (than shown) S is often oscilliatory in na¬ 
ture. In these simulations at least, when the nonlinearity is 
strong enough to induce a shock, it happens after an interval 
of uniformly increasing S. 

An alternative behaviour is shown on fig. where the 
progress of the shock coefficient 5(f) changes with an increas¬ 
ingly group velocity Vg. We can see that once the group ve¬ 
locity becomes significant, the maximum wave gradients as 
indicated by 5 drop quickly, and shocking no longer occurs; 
the steepness also oscillates as time passes. 

Lastly, fig. shows the progress of the shock coefficient 
5(f) as it changes with an increasing dispersion d 2 . We can 
see that as is not unexpected ||7]|9l, the response to ever in¬ 
creasing dispersion tends to be a straightforward monotonic 
reduction in the maximum wave gradient. However, the inter¬ 
play of nonlinearity, group velocity, and dispersion can some- 


FIG. 3: Shock progression in a dispersionless medium, with the 
simulation data used to create from fig. 1(a) of Kinsler Q. (a) Wave 
profile, where the solid, dashed, dash-dotted, and dotted lines are for 
propagation distances closer and closer to the shocking distance, (b) 
Shocking coefficient for the same data, which increases as expected 
and also indicates the moving location of the most shocked part of 
the wave. 



FIG. 4: Shock coefficient 5(f) for a set of x^^'^ media with vary¬ 
ing nonlinearity. The basic simulation is the same as for fig. 1 from 
Genty et al. (6), but with fixed group velocity Vg = 1.00 and disper¬ 
sion d 2 = 0.50. The “missing” 5(f) data at longer times and larger 
nonlinearities occurs because an LDD shock was detected and the 
simulation terminated. 
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FIG. 5: Shock coefficient S{t) for a medium with a varying 
group velocity v^. The basic simulation is the same as for fig. 
above, but with a fixed nonlinearity of = 0.32 and disper¬ 

sion d 2 = 0.50. 



FIG. 6: Shock coefficient S{t) for a x^^'^ medium with a varying 
dispersion d 2 . The basic simulation is the same as for fig. |^above, 
but with a fixed nonlinearity of x^^'^ = 0.32 and group velocity 

Vg = 0.60. For a few smaller values of d 2 it happens that an LDD 
shock was detected and the simulation terminated; the lack of data 
thereafter gives rise to the subsequent trough (with an artificial S = 
—0.5) at larger times. 


times cause shocking even when ordinarily it might not be 
expected. Although not shown here, some simulations with 
similar parameters to those used for fig. [^did detect (LDD) 
shocking after short propagation distances, notably for limited 
ranges of the smaller values of d 2 . 

The main utility of the proposed shocking coefficient S is 
that we no longer need to look as figures such as fig. 1 in ifTSll 
or fig. 15 in 0, and by eye try to estimate how much more 
or less a waveform is shocked compared to some compara¬ 
tor. Instead we can calculate S for those waves and generate 
results like those in figs. 00 or|^ Given these, it is straight¬ 
forward to notice that a given change in parameters causes S 
to increase or decrease by some calculable amount; or that S 
has (also) passed some well-chosen threshold value. It may be 
noted that graphs showing maximum gradient of an evolving 
waveform have been presented in the past while the shocking 
process is being discussed CHi, but while these allow us to 
judge how much the field gradient has steepened, they are less 
useful in informing us of how far they are from a shock (i.e. 
from infinite steepness). 

As a final comment, all the simulations shown above were 
done using the computational model of Tyrrell at al. ini; 
this is a full field approach which propagates pulses in space. 
Spatial propagation methods are extremely common in non¬ 
linear optics, but this is necessarily an approximate approach 
nni [El [Ts] m, and caution may be needed in cases with 
extreme nonlinearity or dispersion. However, the definition 
of S used here is not dependent on how a waveform was ob¬ 
tained, and once could as easily evaluate a temporally propa¬ 
gated S{x) as the S{t) considered here. 


IV. CONCLUSIONS 

In this paper I have proposed a quantitative measure of 
shocking, derived from the gradient of the wave profile, with 
the infinite steepness at shocking moderated down to a finite 
value by using the inverse tangent function. The results in 
sec. [^suggest that that this definition of S will enable us to 
avoid the imprecision that follows a more qualitative judge¬ 
ment based solely on (e.g.) the appearance of a wave profile. 
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